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We introduce a coherence susceptibility method, based on the fact that it signals quantum fluctuations, for 
identifying quantum phase transitions, which are induced by quantum fluctuations. This method requires no 
prior knowledge of order parameter, and there is no need for careful considerations concerning the choice of 
a bipartition of the system. It can identify different types of quantum phase transition points exactly. At finite 
temperatures, where quantum criticality is influenced by thermal fluctuations, our method can pinpoint the 
temperature frame of quantum criticality, which perfectly coincides with recent experiments. 


INTRODUCTION 

Fluctuations trigger phase transitions manifesting them¬ 
selves as the sudden change of the system states in the ther¬ 
modynamic limit [1-3]. At absolute zero temperature ther¬ 
mal fluctuations cease, leaving quantum fluctuations the only 
source for the corresponding transitions, which are known as 
quantum phase transitions (QPTs). Arising from the Heisen¬ 
berg uncertainty principle quantum fluctuations underlie the 
“quantumness” of the many-body systems. 

Different quantum phases of matter had been in general 
characterized by their symmetries under Landau’s symmetry- 
broken theory before the striking discovery of the (fractional) 
quantum Hall effect in the 1980s [4]. Since then exotic quan¬ 
tum phases such as topological ordered phases have emerged 
as one major topic in condensed matter physics and many- 
body physics. In modern times, with the development of 
quantum information sciences a wide variety of characteriza¬ 
tions for quantumness have been proposed for the investiga¬ 
tion of QPTs including topological phase transitions. Entan¬ 
glement was the first and most famous one [5-9]. Quantum 
discord, which measures the quantum correlation between two 
components of the system, complemented entanglement in 
certain situations to detect QPTs [10, 11]. As another indica¬ 
tion of quantumness, i.e., the lack of local convertibility, was 
also useful in characterizing quantum phases [12-15]. Com¬ 
pared with traditional methods based on condensed matter 
physics, which underlie low-energy effective theories of lo¬ 
cal order parameters [16], the quantum information oriented- 
methods are provided with a common advantage in that they 
can study QPTs without any knowledge of order parameters 
in priority. This merit has been validated in particular through 
the studies of exotic quantum phases with topological orders 
[15, 17-21]. Though remarkable success has been achieved, 
some obvious drawbacks exist, e.g., a careful consideration 
of the bipartition of the system is required by definition. For 
instance, it has been shown in Ref. [22] that using the entan¬ 
glement method some OPTs can be detected only by bipartite 
entanglement between certain two sublattices but not by pair¬ 
wise concurrence between two nearest-neighbour sites. 


Another important topic in studying QPTs is to And out the 
influence of thermal fluctuations on the quantum criticality of 
QPTs, since in the real world QPTs have to be observed and 
manipulated at finite temperatures. It is known that certain 
QPTs or quantum criticality persists at non-zero temperatures 
[3, 5, 23-26]. For instance, recently the temperature frame 
of quantum criticality of an Ising model QPT has been mea¬ 
sured experimentally [26]. The quantum criticality has been 
reported to survive up to temperature T ~ 0.4J with J be¬ 
ing the coupling strength of the Ising model [26]. The closest 
theoretical prediction was T ~ 0.5 J [24], which is not so de¬ 
sirable. 

Both aspects above are related to the coherence of the sys¬ 
tem. Indeed, coherence is a better alternative to depict quan¬ 
tumness as it involves no partition of the system. Coherence 
has been widely used in the fields of biological physics [27] 
and quantum open systems [28]. It was mainly characterized 
by the off-diagonal entries of the density operator. The prob¬ 
lem was that the many measures of coherence were merely out 
of physical intuitions. Only recently has a rigorous framework 
been established [29] where the measure based on relative en¬ 
tropy and that on the li matrix norm were shown as legitimate 
measures. Much attention has been paid in this direction since 
then [30-33]. The basic idea in Ref. [29] is to treat coher¬ 
ence under the resource theory, and the central requirement 
is that coherence should not increase under incoherent oper¬ 
ations. The following simplified form of relative entropy has 
been proven as a valid measure of coherence for a given basis: 

C{p) = S{pdiag) - S{p), (1) 

where S{») stands for the von Neumann entropy of • and 
Pdiag is obtained from p by removing all its off-diagonal en¬ 
tries. 

On the other hand, QPTs can be formulated as follows. A 
many-body system with tunable parameter(s) A is described 
by the Hamiltonian H{X). For each given A the ground state 
of the system is labeled as |T'g(A)). At QPT point Ac when 
slightly varying A, the ground state drastically changes. Oc¬ 
curring at absolute zero where thermal fluctuation is com¬ 
pletely frozen, QPT roots in the quantum fluctuation only. 
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which is induced by the change of A. Hence, the non- 
analyticity of the ground state at phase transition point can be 
characterized by the singularity of the coherence susceptibil¬ 
ity, which is defined as s dC{p)ldX. Here, p can be the 
density operator of the whole system or the reduced density 
operator of a subsystem. Throughout this paper the deriva¬ 
tive in the coherence susceptibility is carried out numerically, 
whereas the coherence is calculated analytically because of 
the exact solutions of the corresponding models. 

In this paper, we will establish the relative entropy in Eq.(l) 
as the coherence measure and its susceptibility as a powerful 
tool to explore QPTs. The coherence susceptibility method 
reflects the origin of QPTs, i.e., quantum fluctuations, and 
it requires no prior knowledge of the order parameters nor a 
careful consideration of the bipartition. Its simplicity makes it 
suitable to study QPTs in general, and to investigate the com¬ 
plicated situation of finite-temperature QPTs. We showcase 
the performance of our method through several models. The 
temperature frame of quantum criticality pinpointed by our 
method perfectly matches the experimental result in the Ising 
model [26]. 


RESULTS FOR ZERO TEMPERATURE 

We apply the coherence susceptibility method to two spin 
models with different types of ordinary QPTs and one model 
with topological QPT Since the measure of coherence de¬ 
pends on the basis of writing the density matrix, we fix the 
computational basis throughout the paper, i.e., the product of 
eigenbasis of to write the reduced density matrices of the 
ground states for these models and then substitute p in Eq.(l) 
with them to calculate the coherence. Note that choosing a 
different basis may change the value of coherence (and its 
susceptibility) for each state in the phase diagram; however, 
the singularity stays at the phase transition point. In very rare 
cases, the coherence may always equal zero for all states by 
poorly choosing certain basis, where no useful information 
can be extracted. Such a negative fine-tuned effect can be 
avoided in general by taking another basis. One also has to 
pay attention to the size of the reduced density operator. In 
general the bigger the reduced density operator, the more in¬ 
formation it contains, which means a higher chance to probe 
the QPT but being more difficult to calculate. In practice one 
can start from the one-site density operator to calculate the 
coherence. If the result is trivial, meaning that the coherence 
is always zero and no information can be extracted, one has 
to increase the size of the density operator until a non-trivial 
density operator emerges. 

The transverse field Ising model (TEIM) in one spatial di¬ 
mension is the simplest model with a quantum phase transi¬ 
tion. Its Hamiltonian reads 

= ( 2 ) 

i=l i=l 

where cr®’* are the usual Pauli matrices. We can choose the 
free parameter A := JfB to study the phase diagram. It can be 


solved analytically using Jordan-Wigner transformation [34- 
36]. The one-site reduced density operator for the ground state 
can be recovered from the expectation values of the Pauli ma¬ 
trices as Pi = where and 

= / is the identity matrix. Here (•) denotes the average of 
• over the ground state. The critical point locates at the second 
order phase transition point Ac = 1. Eor 0 < A < 1, the system 
is gapped, and the ground state is in the paramagnetic phase 
with vanishing order parameter {a^)- When A > 1 the system 
is also gapped with double degeneracy in the energy spectrum. 
In the thermodynamic limit the symmetry can be broken spon¬ 
taneously, and the (degenerate) ground states are in the ferro¬ 
magnetic phase with non-vanishing order parameter (cr^). The 
difference between the two branches of the degenerate ground 
states is the opposite sign of the order parameter, and there¬ 
fore they are labeled as |0^) and |0“), respectively. In this 
phase the system is also possible to stay in an equal mixture 
of |0^) and |0~), i.e., ptg = (|0'*') (0'*'| + |0“) (0“|)/2, which 
has zero order parameters. This mixture state is known as the 
“thermal ground state.” In this simple model we can investi¬ 
gate the one-site reduced density operator for the system in 
the thermodynamic limit to detect the QPT. We study the co¬ 
herence for both the symmetry broken ground state and the 
thermal ground state. 

It can be deduced from Ref.[8] that for |0^) 




0 ,A<1, 

(i-a-2)1/8 ,A>1, 


(3) 


1 r'^ 

(a“) = - / # 

TT Jo 




1 -H A cos (j) 

+ A^ + 2Acos(/i 


(4) 


and (cr^) = 0. Eor the thermal ground state the only difference 
is that (tr^) is always zero. Inserting these results to Eq.(l) 
one can obtain the coherence for TEIM in Eig. 1. The main 
plot shows that the quantum phase transition is detected by the 
singularity of the coherence for the symmetry-broken ground 
state, and the inset shows that it is identified by the divergence 
of the coherence susceptibility for the thermal ground state. 

A complementary type of quantum phase transition is that 
in one phase the system is gapped whereas in the whole region 
of the other phase the system is critical. A spin half XX 
model is one such example with Hamiltonian 




1 ^ 

= --E 

2&1 


M+l 




N 

•AE' 

i=l 


(5) 


where the dimensionless parameter A denotes the strength of 
the external magnetic field in units of the interaction energy. 
The phase diagram is symmetric with respect to A [37-39], 
therefore we consider only positive A. This model can be 
solved analytically [37, 39]. At A = 1 the system under¬ 
goes a first order quantum phase transition. When A > 1 
the ground state is polarized up for all spins resulting in van¬ 
ishing coherence. In the region 0 < A < 1 the system is 
critical and the one-site reduced density operator is diagonal. 
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FIG. 1. Ising model ground state. Main plot: Coherence of one-site 
reduced density matrix for the symmetry-hroken ground state (dash- 
dot line) and the thermal ground state (solid line). For the symmetry 
broken state the singularity occurs exactly at the quantum phase tran¬ 
sition point. Inset: coherence susceptibility for the thermal ground 
state. This is obtained by numerically carrying out the derivative of 
the coherence curve in the main plot with respect to A, therefore the 
actual divergence of x at the phase transition point manifests itself as 
a big negative value in the numerical result due to the finite step. 

so its coherence always vanishes according to Eq.(l). The 
one-site reduced density operator is trivial in terms of co¬ 
herence since no information can be extracted. We can em¬ 
ploy the coherence of two adjacent spins to detect the quan¬ 
tum phase transition point. Its reduced density operator reads 
P 2 = Za=o E|=o ® ® cr^/4. It has been shown 

in Refs. [37, 39] that the non-trivial non-vanishing coeffi¬ 
cients are (cr^cr^) = (l - 2arccos(A)/7r)^ - 4(1 - A^)/7r^, 
(tJ^cr®) = -2sin(arccos(A))/7r and (cr^) = l-2arccos(A)/7r. 
Substituting these in Eq.(l) one obtains the coherence in Eig. 
2, where the singularity at A = 1 perfectly detects the phase 
transition. 

There exists another kind of exotic quantum phase transi¬ 
tions, namely, topological quantum phase transitions, without 
any local order parameter and therefore they cannot be de¬ 
scribed by the Landau paradigm [40] . We will show through 
the Kitaev honeycomb model that the topological phase tran¬ 
sition can be captured by the coherence susceptibility method. 
The Kitaev honeycomb model on a hexagonal lattice with 
direction-dependent interactions between adjacent lattice sites 
is an analytically solvable model with topological quantum 
phase transition [41]. Its Hamiltonian reads 

hk = - z E (6) 

oc={x,y,z} (i,j)eQ;-Zm/cs 

where the parameters Ja represent the interaction energy be¬ 
tween two adjacent spins along the a direction with the a link. 
We fix Ja; + Jj, + = 1 as the energy unit for this model. 

The phase diagram is shown in the left inset of Eig. 3. In 



FIG. 2. The coherence of two adjacent spins for the XX model in the 
thermodynamic limit. 

the shaded areas the system is gapped with Abelian excita¬ 
tion, and in area B it is gapless with non-Abelian excitation. 
We take the path Jy = = (1 - Ja;)/2 to study the topo¬ 

logical phase transition, which occurs at = 0.5. The an¬ 
alytical solution of the ground state is applied to derive the 
reduced density matrix of two adjacent spins with x link in 
the thermodynamic limit. The only non-vanishing non-trivial 
coefficient is {a^a^) = + A^dwycitLi^/dTr^ with 

e = Jx + Jy cos LOy + Jz cos iOz and A = Jy sin ujy + J^ sin uj^ 
(see the Appendix for details). For this model, the coherence 
susceptibility is defined as % = dCIdJx- Figure 3 depicts the 
two-site coherence and its susceptibility for the Kitaev model. 
The topological quantum phase transition at J^; = 0.5 is well 
captured as a singular point. 


RESULTS FOR FINITE TEMPERATURES 

In real experiments, finite temperature is unavoidable due 
to the third law of thermodynamics, hence the influence of 
thermal fluctuation has to be considered. At finite tempera¬ 
tures, quantum phase transitions are washed out by the ther¬ 
mal fluctuations, nevertheless certain phase transitions still 
persist in 2D systems [3]. For ID systems, e.g., the TFIM, 
the singularity of the free energy disappears at finite tempera¬ 
tures [23] however, the quantum phase transition has so pro¬ 
found an influence that the criticality is detectable even at 
finite temperatures [5, 24, 25]. In other ID examples, the 
phase transition points at zero temperature were reported to 
be signalled by the quantum correlations at finite tempera¬ 
tures [11]. In this context the system to study is in thermal 
equilibrium with a reservoir at temperature T. Thus, it is de¬ 
scribed by the Gibbs thermal state as p = exp{-HlkBT)lZ, 
where Z = Ti{exp{-H/kBT)} is the partition function, kB 
is the Boltzmann constant, and sets the energy scale of 
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FIG. 3. Kitaev model ground state. The left inset is the phase dia¬ 
gram of the Kitaev model. We calculate the ground state along the 
dashed line. The right inset shows the obtained coherence for two 
sites with an x link. The main figure plots the coherence suscep¬ 
tibility. The singular point occurs exactly at the topological phase 
transition point Jj; = 0.5. 


the thermal fluctuations. Using TFIM [Eq.(2)] at ksT rang¬ 
ing from 0 to 1 as an example, we showcase here that the 
coherence susceptibility is a powerful tool to study the finite 
temperature phase diagram. In order to compare with the pre¬ 
vious results in this literature, we rewrite the Hamiltonian as 

N N 

Hi = - XI cf <8 - A XI erf- The finite temperature phase 

i=l i=l 

diagram was predicted in Ref. [25] [also see the left inset of 
Fig. 4 where the two dotted crossover lines separate it into 
three distinct regions; renormalized classical (RC), quantum 
critical (QC), and quantum disordered (QD) region]. On the 
crossover lines, k^T equals the energy gap, i.e., 2|A-Ac|. Re¬ 
cently this diagram has been confirmed in experiments [26]. 
Moreover, the experiment showed that the criticality could 
survive up to ksT rs 0.4. In the following we show that all 
the above complex properties can be detected by the simple 
coherence susceptibility method. 

Again, we take the analytical form of the one-site reduced 
density operator for the thermal state with 

K).-i (7) 

TT w/o 

where uj^ = ^Jl + + 2Acos (p, and (u^) = (u^) = 0 for the 

thermal state [42]. Fig. 4 demonstrates its coherence sus¬ 
ceptibility. One can And that with growing temperature the 
singular behavior deforms from the divergence to maximum, 
which shows that the quantum phase transition is smoothed 
out by the increasing thermal fluctuations. The right inset ac¬ 
counts for the location of the maximum x, which is labeled 
as Xm{T), at different temperatures. When kgT is below 0.4 
the Xm{T) is linear with T, which indicates the criticality of 
the model. The linear scaling changes at ksT = 0.4, which 



A, 



FIG. 4. Ising model thermal states at different temperatures. Main 
plot: coherence susceptibility at different temperatures. At finite 
temperatures the quantum phase transition is smoothed out due to 
thermal fluctuation. As a result the singularity deforms into maxi¬ 
mum points. Left inset: Phase diagram at finite temperatures. Right 
inset: the locations of the maximum point at different temperatures. 
Here fcsT is in unit of the interaction energy in the system Hamilto¬ 
nian. 


coincides with the temperature frame in experiment [26]. If 
one compares the Xm{T) curve with the phase diagram in the 
left inset, it is clear that when fc^T < 0.4, the Xm{T) are lo¬ 
cated on the second crossover line kgT = 2(A - Ac), so that it 
not only detects the temperature range of criticality, but also 
pinpoints the crossover boundary between the QC and QD re¬ 
gions. From the experimental point of view, this feature is 
useful to extrapolate the quantum phase transition point (oc¬ 
curring at T = 0) with data from real experiment (at finite 
temperatures). 

To compare, in Ref. [11] quantum discord was reported 
to be able to exactly detect the quantum critical point of the 
XXZ model even at finite temperatures; however, we found 
that quantum discord is unable to detect the temperature frame 
for the quantum criticality of Ising model. Since quantum dis¬ 
cord is invariant under local unitary transformations, we can 
change the basis by exchanging and on every site. Then 
the reduced density matrix of the two nearest neighbor sites 
can be written as p 2 = [/ 8 / + {a^a^) + {a^a^) 8 

a^ + {a^cr^) a^^a^+{a^) which is of the 

“X” type so that the analytical solution for quantum discord in 
Ref [43] can be applied, (cr^) is already known from Eq.(7), 
and the others can be easily derived from Ref. [42] as 


where 


{ayay) = G,, 

(aV^) = G_, 

1 ,, (A-H cos(/)) tanh(-w^/fc_BT) 

— / dej) - 

TT Jo 

±— / d(j) -, 

TT Jo 


( 8 ) 

(9) 

( 10 ) 


( 11 ) 


with uj^ = y/l + X^ + 2AcoS(;i. 
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FIG. 5. Quantum discord for Ising model thermal states at different 
temperatures. From top to bottom the temperature increases from 
fcsT = 0 to 1 in unit of the interaction energy in the system Hamil¬ 
tonian. 



FIG. 6. The A location of the maximal quantum discord with re¬ 
spect to temperatures for the Ising model. We show it only up to 
ksT = 0.7 because the real maximal quantum discord for higher 
temperatures is out of the calculation scope of A. Here ksT is in 
unit of the interaction energy in the system Hamiltonian. 


DISCUSSION 


We have introduced a coherence susceptibility method to 
detect general quantum phase transitions. Different from 
the entanglement method and quantum discord method [11], 
which succeed in detecting QPTs essentially by selecting var¬ 
ious correlation functions in some educated ways, the under¬ 
lying idea of the present paper is that coherence represents 
the quantumness of a system, and its susceptibility reflects the 
quantum fluctuation, which triggers the quantum phase transi¬ 
tion. As a result at quantum phase transition points the coher¬ 
ence susceptibility should show singularities. The advantage 
of this method over other existing methods is that it requires 
no prior knowledge of the order parameters of the system and 
no careful consideration of the bipartition is needed. In this 
sense, it is akin to the fidelity susceptibility method but with 
a different mechanism as its foundation. We have applied this 
method to various quantum systems with continued, first or¬ 
der, and topological phase transitions. In all cases, our method 
can detect the phase transition points exactly. 

The effect of finite temperature on the coherence suscepti¬ 
bility method has also been discussed. Our method not only 
can extrapolate the quantum phase transition point with finite 
temperature data, but also pinpoint the crossover boundary be¬ 
tween quantum critical and disordered regions and estimate 
the temperature frame of quantum criticality. The tempera¬ 
ture frame predicted by our method matches the experimen¬ 
tal conclusion perfectly, which has filled the gap between the 
previous theoretical result [24] and the recent experiment. It 
is an interesting open question as to whether the crossover 
boundary that separates the quantum critical region and the 
renormalized classical region can be detected through the co¬ 
herence susceptibility of bigger constituent systems. One step 
further would be to generalize the coherence susceptibility 
method from identifying quantum phase transition points to 
the quantum-classical crossover[44, 45]. It is also interesting 
to investigate the dynamics of closed or open systems at both 
zero and finite temperatures using the coherence susceptibil¬ 
ity. 
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APPENDIX: THE EXACT SOLUTION OF THE KITAEV 
HONEYCOMB MODEL 

In this appendix we show how to derive the ground state of 
the Kitaev model and calculate the correlation function of two 
sites with an x link using the original method in Ref. [41]. 

We first introduce the Majorana operators. At each site, we 
define four Majorana operators c“, with a = 0,x,y,z, satis¬ 
fying {c°‘y = c“, {c“,c^} = 2Sai3, and = 1. We 

denote c° by Cj and represent the Pauli operators by the Ma¬ 
jorana operators as 

(a = a:,2/,z) (12) 

Since the values of a in the Kitaev model are determined by 
the site index j and k, we can rewrite the Hamiltonian as 

^ Jaj^k'^j,kCjCk iUj,k = iCjC).) (13) 
^ j,k 

Kitaev has showed that = 1, [uj^k, H]=Q, and com¬ 
mute with each other. We take Uj^k = 1 for all links because 
this vortex-free configuration has the lowest energy [46]. We 
follow the convention that j is the sublattice presented by the 
empty circles and k is the other sublattice presented by the full 
circles in Fig. 7. 



FIG. 7. A graphic representation of the Kitaev model. There are 
two sublattices ( empty and full circles). Each unit cell (marked by 
elliptic circle) contains one site of each kind. Three types of bonds 
are labeled by “x-links,” “y-links,” and “z-links”. We choose the 
coordinate axes in ni and n 2 directions. 

Then, we choose the unit cell containing one empty circle 
and one full circle along the x link, shown in the Fig. 7. We 
rewrite the site index j as (s, A), where s refers to the location 
of the unit cell, and A describes the two different kinds of 
sublattices(empty circles take the value 1, full circles take the 
value 2). Thus, the Hamiltonian becomes 

i 

“ o 'fs,A,t,/2Cs,AQ,/x- (14) 

S,A,t,/L 

Actually, Js,A,t,/i is determined by three indexes A, p, and t-s. 
Then, using the Fourier transformation, the Hamiltonian takes 


the form 

2 __ 

* X! S (15) 

q 

with 

t 

Oq.A = E e”*'*’''’Cs,A, (16) 

where a ,,a satisfies a_q,A = = 0 ’ {ap,A>«q.M} “ 

6pq6\^fj,, and other anticommutators are all equal to zero. 

After simple calculations we obtain that A.i(q) 
</ 2 , 2 (q) = 0. We choose two directions ni = V3/2ey - e^/2, 
and n 2 = \/3/2ej, + 62 / 2 , shown in Fig. 7. Then, Ji, 2 (q) = 
Jx + + J 2 e*‘^■"^ J2,i(q) = -Jl^{q). Let /(q) = 

e(q) + iA(q), and choose to be in the direction of ni, and 
^ to be in the direction of n 2 . Then we have 

e(q) = Jx + 4 cos qy + Jz cos qz, 

A{q) = Jysmqy +JzSinqz, (17) 

where qy and q^ take values qy,qz = 27rn/L, n = -{L - 
l)/2 ,-,(L-l)/2 . 

Next, we use the Bogoliubov transformation to diagonalize 
the Hamiltonian, 

61q,l — tXqtlq 1 + 'UqUq, 2 ; 

C t * t . ’I- t 

q,l = ^^q<l +^^q«q.2. 

61q,2 “ t^qtTq,! “ t/qUq^ 2 j 

Cq .2 ( 1 ^) 

where, Uq = l/^/2, Vq = ifq/{^/2\fq\), u_q = -Uq, and 
the new operators satisfy {C'q.AjCp^x } = ^pq<5A,;., Cl^ = 0 , 

C_q,i = -2uy^Cl^, and C^Cq,! = 1 - Cl, 2 C_q, 2 . Then, 
the Hamiltonian reads 

7T=^|/q|(l-2CEC'q,2). (19) 

q 

The normalized ground state is 

|G') = ncElO), (20) 

q 

with Cq,2|0) = 0. The energy gap is 2minq{|/q|}. 

CORRELATION FUNCTIONS AND REDUCED DENSITY 
MATRIX OF THE KITAEV MODEL 

The two-site reduced density matrix is the joint state 
of two spins at sites i and j and takes the form pi j = 
Z'f’ where cr° is the identity, and = 

need 16 correlation functions to construct the 
two-site density matrix. We calculate the correlation functions 
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of two nearest lattices linked by an x bond, and it is obvious 
that 1 ^^ 2 ) - 1- For rhe rest of the correlation functions, 
only (tTp iCTf 2 ) is nonzero. Thus, the two-site reduced density 
matrix along the x link has the form. 


P2 = 


/ 1 
0 
0 

\ {a^a^) 


0 

1 


(cr^CT^) 


{a^a^ 

0 

0 

1 


)\ 


( 21 ) 


Then, we calculate the correlation function (<J^iO'r 2 ) using 
the analytic result of the model: 


« 1 < 2 ) = E ^■'■aq.lV, 2 ) 


q q 




q.q 


By using the relation 


(Uq,lUq',2) “ ((UqCq,! + UqC'q^2 ) (u^/Uq'Cq' 2 )) 
= -UqMq'(Cq^Cq' 2) 


= -6 -A 

“2 l/ql’ 

we obtain the correlation function 


/q 


1 V- A + 


«1<2) = ^2 E = 2L2 S 

- — ^JL 


( 22 ) 


where T;q = |/q| = + A^. 

In the thermodynamic limit, we use the continuous u!y(^z) 
replacing qy(z), then the congelation function has the form 


■)■—[' r 

' 47r^J-, J- 




dujyduiz, ( 23 ) 


with e = Jx + Jy cos ujy + Jz cos ujz and A = Jy sin ujy + Jz sin^. 
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